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1. Introduction 

The QCD transition at non-vanishing temperatures (T) plays an important role in the physics 
of the early Universe and of heavy ion collisions (most recently at RHIC at BNL; LHC at CERN 
and FAIR at GSI will be the next generation of machines). The main goal of the present summary 
is to present results of the Budapest- Wuppertal group on the QCD transition at vanishing chemical 
potential (jli=0), which is of direct relevance for the early universe (/I is negligible there) and for 
present and future heavy ion collisions (at RHIC £i<40 MeV, which is far less than the typical 
hadronic scale). Since these results were obtained at four different sets of lattice spacings and a 
careful continuum extrapolation was performed, we consider them as full results. In addition, a 
new technique will be introduced, which closes the gap between lattice QCD and the perturbative 
approach for bulk thermodynamical quantities. 

The standard picture for the QCD phase diagram on 
the light quark mass («w) versus strange quark mass (m s ) 
plane is shown by Figure [I]. It contains two regions at 
small and at large quark masses, for which the T > 
QCD transition is of first order. Between them one finds 
a cross-over region, for which the T > QCD transition 
is an analytic one. The first order transition regions and 
the cross-over region are separated by lines, which corre- 
spond to second order phase transitions. 

When we analyze the nature and/or the absolute scale 
of the T > QCD transition for the physically relevant 
case two ingredients are quite important. 

First of all, one should use physical quark masses. 
As Figure [I] shows the nature of the transition depends on 
the quark mass, thus for small or large quark masses it 
is a first order phase transition, whereas for intermediate 
quark masses it is an analytic cross over. Though in the 
chirally broken phase chiral perturbation theory provides 
a controlled technique to gain information for the quark mass dependence, it can not be applied for 
the T > QCD transition (which deals with the restoration of the chiral symmetry). In principle, 
the behaviour of different quantities in the critical region (in the vicinity of the second order phase 
transition line) might give some guidance. However, a priori it is not known how large this region 
is. Thus, the only consistent way to eliminate uncertainties related to non-physical quark masses is 
to use physical quark masses (which is, of course, quite CPU demanding). 

Secondly, the nature of the T > QCD transition is known to suffer from discretization er- 
rors [|], ^, ]3|]. The three flavour theory with standard action on N t =4 lattices predicts a critical 
pseudoscalar mass of about 300 MeV. This point separates the first order and cross-over regions of 
Figure 1. If we took another discretization, with another discretization error, e.g. the p4 action and 
N t =4, the critical pseudoscalar mass turns out to be around 70 MeV (similar effect is observed if one 
used stout smearing improvement and/or N t =6). Since the physical pseudoscalar mass (135 MeV) 
is just between these two values, the discretization errors in the first case would lead to a first order 
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Figure 1: The phase diagram of QCD 
on the hypothetical light quark mass ver- 
sus sttange quark mass plane. Thick lines 
correspond to second order phase transi- 
tions, the purple regions represent first or- 
der phase transitions and the yellow re- 
gion represents an analytic cross-over. 
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Figure 2: The ratio f(cont)lf{N t ) as a function of l/N} (left panel). f(cont) is the continuum extrap- 
olated free energy of the staggered fermionic gas in the non-interactive, infinitely high temperature limit. 
f{Nt ) is the value obtained on a lattice with N t temporal extension. The black line shows our choice (stout 
improvement, only next-neighbour terms in the action), whereas the red and blue lines represent the Naik 
and p4 actions, respectively. Masses and taste symmetry violation for different approaches in the literature 
(right panel). The smallest, physical quark mass and the smallest taste symmetry violation was reached by 
our works (black dot, ||, |J]). Somewhat larger taste symmetry violation and about three times larger quark 
masses were reached by the MILC analysis on QCD thermodynamics (red dot, JTc|]). Even larger taste 
symmetry violation and about four times the physical quark masses are the characteristics of the Bielefeld- 
Brookhaven-Columbia-Riken result on T c (blue dot, |Q). 

transition, whereas in the second case to a cross-over. The only way to resolve this inconclusive 
situation is to carry out a careful continuum limit analysis. 

Since the nature of the transition influences the absolute scale (T c ) of the transition -its value, 
mass dependence, uniqueness etc.- the above comments are valid for the determination of T c , too. 

Therefore, we have to answer the question: what happens for physical quark masses, in the 
continuum, at what T c ? In order to get a reliable answer we used physical quark masses on N t =4,6,& 
and 10 lattices, which correspond to approximately 0.3, 0.2, 0.15 and 0.12 fm lattice spacings, 
respectively. 

In the presentation [Q] of the results of the Bielefeld-Brookhaven-Columbia-Riken Collabora- 
tion published results [Q] from N r =4 and 6 were shown (and some unpublished figures for N t =8, 
which were obtained within the HotQCD Collaboration). Since the CPU requirements for thermo- 
dynamics increase as as our N t =\0 simulations need about 50 times more CPU than N t =6. Do 
we have 50 times more resources for QCD thermodynamics than our competitors? Of course not 
(it is almost the other way around). Instead, reaching N t =\0 is a fine balance. It is partly related to 
the choice of our action (which will be discussed in the next section), partly to the arrangements of 
the financial resources. For instance, as N t increases, one needs more and more statistics. Thus the 
thermalization can be done only once on a relatively expensive, scalable machine, such as Blue- 
Gene/L, whereas a large fraction of the non-vanishing T simulations can be done on more cost 
effective devices such as personal computer [Q] graphics cards. A 2 years old model can accommo- 
date N t =6 lattices, on a one-and-a-half year old model you can put A^=8 lattices and the one year 
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Figure 3: The volume dependence of the susceptibility peaks for pure SU(3) gauge theory (Polyakov-loop 
susceptibility, left panel) and for full QCD (chiral susceptibility on N t =4- and 6 lattices, middle and right 
panels, respectively). 

old model can work with quite large N t =lO lattices. It costs a few hundred dollars and can provide 
upto 30-60 Gflops sustained QCD performance. They are not easy to code, adding two numbers 
needs 3 pages, but recently more efficient programming environments were introduced. Clearly, 
this type of hardware provides a very advantageous price-performance ratio for lattice QCD. 

2. The choice of the action 

The first step is to choose an action, which respects all the needs of a thermodynamic analysis. 
T=0 simulations are needed to set the scale and for renormalization. T>0 simulations are needed 
to map the behaviour of the system at non-vanishing temperatures. The action should maintain a 
balance between these two needs, leading to approximately the same uncertainties for both sectors 
(otherwise a large fraction of the CPU-power is used just for "over-killing" one of the two sectors). 
We used Symanzik improved gauge and stout improved 2+1 flavour staggered fermions [f7|] (due 
to the stout improvement we have only next-neighbour terms in the fermionic part of the action). 
The simulations were done along the line of constant physics. The parameters were tuned with 
a quite high precision, thus at all lattice spacings the nin/fk and m/f/m^ ratios were set to their 
experimental values with an accuracy better than 2%. 

The choice of the action has advantages and disadvantages. As we will see the advantages are 
probably more important than the disadvantages. The left panel of Figure |2] shows the continuum 
free energy devided by its value at a given N t . A related plot is usually shown by the Bielefeld- 
Brookhaven-Columbia-Riken collaboration as a function of N t . Since in staggered QCD most 
lattice corrections scale with a 2 , which is proportional to l/N 2 , it is instructive to show this ratio as a 
function of l/N 2 , for our action, for Naik and for p4. Extrapolations from 4 and 6 always overshoot 
or undershoot. Clearly, the Naik and p4 actions reach the continuum value much faster than our 
choice, but the a 2 scaling appears quite early even for actions with next-neighbour interactions. 
Extrapolations from N t and N t +2 with our action are approximately as good as N t with the p4 
action (which was taylored to be optimal for this quantity, namely for the free energy at infinitely 
large temperatures). In practice, it means that our choice with _/V,=8,10 gives approximately 2% 
error for the free energy. In a balanced analysis you do not need more, because the corresponding 
lattice spacings 0. 15 and 0. 12 fm are most probably not fine enough to set the scale unambiguously 
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Figure 4: Continuum extrapolated susceptibilities T 4 /(m 2 Ax) as a function of 1/(T^V). For true phase 
transitions the infinite volume extrapolation should be consistent with zero, whereas for an analytic crossover 
the infinite volume extrapolation gives a non-vanishing value. The continuum-extrapolated susceptibilities 
show no phase-transition-like volume dependence, though the volume changes by a factor of five. The V^°° 
extrapolated value is 22(2) which is 11a away from zero. For illustration, we fit the expected asymptotic 
behaviour for first-order and 0(4) (second order) phase transitions shown by dotted and dashed lines, which 
results in chance probabilities of 10~ 19 (7 x 10~ 13 ), respectively. 

with the same accuracy. (E.g. the asqtad action at N t ^lO, which corresponds about a=0. 12 fm 
lattice spacing, still has «10% scale difference between r\ & fx-) Since the p4 action is almost 20 
times more expensive than the stout action, it is not worth to pay this price and improve one part of 
the calculation, which hinders you to reach a reasonable accuracy in another part of the calculation. 

This is the balance one should remember in thermodynamics. Indeed, taste symmetry violation 
should be suppressed for many reasons (setting the scale at T = 0, restoring chiral symmetry at 
T > etc.) As it was argued above it is more important to improve on this sector of the calculation 
than on the infinitely high temperature behaviour. The left panel of Figure || shows the splitting 
between the Goldstone and the first non-Golds tone pions for the Bielefeld-Brookhaven- Columbia- 
Riken Collaboration (which is far beyond their kaon mass), for MILC and for us. Our small splitting 
is partly related to the stout improvement and partly to the cost issues, since smaller lattice spacings 
could be used, which resulted in smaller splitting. In addition, the cost issues allowed us to use in 
our finite T simulations physical quark masses instead of much larger masses. 

3. The nature of the QCD transition 

The next topic to be discussed is the nature of the QCD transition. Physical quark masses were 
used and a continuum extrapolation was carried out by using four different lattice spacings. The 
details of the calculations can be found in [g]. In order to determine the nature of the transition one 
should apply finite size scaling techniques for the chiral susceptibility % = (T /V) ■ (d 2 \ogZ / dm 2 ud ) . 
This quantity shows a pronounced peak as a function of the temperature. For a first order phase 
transition, such as in the pure gauge theory, the peak of the analogous Polyakov susceptibility 
gets more and more singular as we increase the volume (V). The width scales with 1/V the height 
scales with volume (see left panel of Figure ||). A second order transition shows a similar singular 
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behaviour with critical indices. For an analytic transition (what we call a cross-over) the peak 
width and height saturates to a constant value. That is what we observe in full QCD on N t =4 and 6 
lattices (middle and right panels of Figure ||). We see an order of magnitude difference between the 
volumes, but a volume independent scaling. It is a clear indication for a cross-over. These results 
were obtained with physical quark masses for two sets of lattice spacings. Note, however, that for 
a final conclusion the important question remains: do we get the same volume independent scaling 
in the continuum; or we have the unlucky case what we had in the Introduction for 3 flavour QCD 
(namely the discretization errors changed the nature of the transition for the physical pseudoscalar 
mass case)? 

We carried out a finite size scaling analyses 
with the continuum extrapolated height of the renor- 
malized susceptibility. The renormalization of the 
chiral susceptibility can be done by taking the sec- 
ond derivative of the free energy density (/) with 
respect to the renormalized mass (m r ). We apply 
the usual definition: f/T 4 =-N? 
■[logZ(N,,N t )/(N t N?) - logZ(N s0 ,N t0 )/(N t0 N^)}. 

This quantity has a correct continuum limit. The 
650 700 750 M J 

t[k] subtraction term is obtained at T=0, for which sim- 

ulations are carried out on lattices with N s q, N t o 
Figure 5: The water-vapor phase diagram. spatial and temporal extensions (otherwise at the 

same parameters of the action). The bare light quark 
mass (mu£) is related to m r by the mass renormalization constant m r =Z m -m u d. Note that Z m falls out 
of the combination m 2 d 2 /dm 2 =m 2 d d 2 /dm 2 u/ . Thus, m 2 d \x(N s ,N t ) — %{Nsa,Nto)\ also has a contin- 
uum limit (for its maximum values for different N t , and in the continuum limit we use the shorthand 
notation m 2 Ax)- 

In order to carry out the finite volume scaling in the continuum limit we took three different 
physical volumes (see Figure ||). The inverses of the volumes are shown in units of T c . For these 
3 physical volumes we calculated the dimensionless combination T 4 /m 2 Ax at 4 different lattice 
spacings: 0.3fm was always off, otherwise the continuum extrapolations could be carried out, 
which are shown on Figure |j. Our result is consistent with an approximately constant behaviour, 
despite the fact that we had a factor of 5 difference in the volume. The chance probabilities, that 
statistical fluctuations changed the dominant behaviour of the volume dependence are negligible. 
As a conclusion we can say that the staggered QCD transition at mu=0 is a cross-over. 




4. The transition temperature 

An analytic cross-over, like the QCD transition has no unique T c . A particularly nice example 
for that is the water-vapor transition (c.f. Figure ||). Up to about 650 K the transition is a first 
order one, which ends at a second order critical point. For a first or second order phase transition 
the different observables (such as density or heat capacity) have their singularity (a jump or an 
infinitly high peak) at the same pressure. However, at even higher temperatures the transition is 
an analytic cross-over, for which the most singular points are different. The blue curve shows the 
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Figure 6: Temperature dependence of the renormalized chiral susceptibility (i^AXyy/T 4 ), the strange 
quark number susceptibility (Xs/T 2 ) and the renormalized Polyakov-loop (Pr) in the transition region. The 
different symbols show the results for N t =4,6,8 and 10 lattice spacings (filled and empty boxes for N t = 4 
and 6, filled and open circles for N t =8 and 10). The vertical bands indicate the corresponding transition 
temperatures and their uncertainties coming from the T^0 analyses. This error is given by the number in the 
first parenthesis, whereas the error of the overall scale determination is indicated by the number in the second 
parenthesis. The orange bands show our continuum limit estimates for the three renormalized quantities as 
a function of the temperature with their uncertainties. 
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a 2 [fm 2 ] 

Figure 7: Difference between the the T c values obtained by 
the Polyakov loop and by the chiral condensate as a function 
of a 2 . 



peak of the heat capacity and the red one the inflection point of the density. Clearly, these transition 
temperatures are different, which is a characteristic feature of an analytic transition (cross-over). 

In QCD we will study the chiral and 
the quark number susceptibilities and the 
Polyakov loop. Usually they give differ- 
ent T c values, but there is nothing wrong 
with it. As it was illustrated by the water- 
vapor transition it is a physical ambigu- 
ity, related to the analytic behaviour of 
the transition. There is another, non-phy- 
sical, ambiguity. If we used different ob- 
servables (such as the string tension, vq, 
the rho mass or the kaon decay constant), 
particularly at large lattice spacings we 
obtain different overall scales. They lead 
to different T c values. This ambiguity disappeares in the continuum limit. According to our ex- 
periences, at finite lattice spacings, the best choice is the kaon decay constant It is known 
experimentally (in contrast to the string tension or ro), thus no intermediate calculation with un- 
known systematics is involved. Furthermore it can be measured on the lattice quite precisely. 

Figure ^ shows the results for the chiral susceptibility, for the quark number susceptibility and 
for the Polyakov loop. Red, blue, green, and purple indicate Af,=4,6,8 and 10 lattices. N t =A is 
always off, the rest scales nicely. The shaded regions indicate the continuum estimates. There is 
a surprising several sigma effect. The remnant of the chiral transition happens at a quite different 
temperature than that of the deconfining transition. It is quite a robust statement, since the Polyakov 
transition region is quite off the ^-peak, and the ^-peak is quite far from the inflection point of the 
Polyakov loop. This quite large differece is also related to the fact that the transition is fairly broad. 
The widths are around 30-40 MeV. 

Due to the broadness of the transition the normalization prescription changes T c , too. It is 
easy to imagine why, just multiply a Gaussian by x 2 and the peak is shifted. That means using 
X/T 2 gives about 10 MeV higher T c than our definition, for which a T 4 normalization was applied. 
(Note, that for the unrenormalized jaT 2 normalization is natural, whereas for the renormalized % 
the natural normalization is done by T 4 . This kind of naturalness manifests itself as possibly small 
errors of the observable.) 

Figure shows the difference between the T c values obtained by the Polyakov loop and by 
the chiral condensate as a function of the lattice spacing squared. The blue band indicates the 
difference for the chiral susceptibility peak position for the T 2 and T 4 normalization. Thus using 
the T 2 normalization no difference can be seen for _/V,=4 and 6, a slight difference is observed for 
^,=8 and a reliable continuum extrapolation needs N t =6,S and 10. 



Our result on T c and that of the MILC Collaboration (7/ c =169(12)(4) MeV [[10|]) are consistent 
within the (quite sizable) errorbars. 

However, our result contradicts the recent Bielefeld-Brookhaven-Columbia-Riken result [0], 
which obtained 192(7)(4) MeV from both the chiral susceptibility and Polyakov loop. This value 
is about 40 MeV larger than our result for the chiral susceptibility (for the Polyakov loop suscep- 



8 



QCD Thermodynamics 



renormalization 
I h 



scale setting [10%] 



overall error 

H 1 



OMeV 



lOMeV 



20 MeV 



30 MeV 



40 MeV 



Figure 8: Possible contributions to the 40 MeV difference between the results of Refs. [M] and 
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Figure 9: Continuum extrapolations based on N t =4 and 6 (left panel: inconsistent continuum limit) and 
using N t =6,8 and 10 (right panel: consistent continuum limit). 



tibility the results agree within about 1.5 a). What are the differences between their analyses and 
ours, and how do they contribute to the 40 MeV discrepancy? The most important contributions 
to the discrepancy are shown by Figure |8[ The first difference is, that in [||] no renormalization 
was carried out, instead they used the unrenormalized quantity %/T 2 . Due to the broadness of the 
distribution this observable leads to about 10 MeV larger T c than our definition. The overall errors 
can be responsible for another 10 MeV. The origin of the remaining 20 MeV is somewhat more 
complicated. One possible explanation can be summarized as follows. In Ref. [||] only Af,=4 and 
6 were used, which correspond to lattice spacings a=0.3 and 0.2 fm, or a _1 =700MeV and lGeV 
These lattices are quite coarse and it seems to be obvious, that no unambiguous scale can be de- 
termined for these lattice spacings. The overall scale in Ref. was set by ro and no cross-check 
was done by any other quantity independent of the static potential (e.g. /#). This choice might lead 
to an ambiguity for the transition temperature, which is illustrated for our data on Figure ^. Using 
only Af,=4 and 6 the continuum extrapolated transition temperatures are quite different if one took 
r or fx to determine the overall scale. This inconsistency indicates, that these lattice spacing are 
not yet in the scaling region (similar ambiguity is obtained by using the p4 action of Having 
^,=4,6,8 and 10 results this ambiguity disappears (as usual N t =4 is off), these lattice spacings are 
already in the scaling region (at least within our accuracy). This phenomenon is not surprising 
at all. As it was already mentioned e.g. the asqtad action at A^IO (which corresponds to about 
a=0.12 fm lattice spacing) has «10% scale difference predicted by r\ or fa- 
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The ambiguity related to the inconsistent continuum limit is clearly unphysical, and it is re- 
solved as we take smaller and smaller lattice spacings (c.f. Figure ||). The differences between 
the T c values for different observables are physical, it is a consequence of the cross-over nature of 
the QCD transition. There is another phenomenon, namely the volume dependence of T c , which is 



also physical. Recently A. Bazavov and B. Berg studied [ |11[ ] the volume dependence of T c in the 
pure SU(3) theory (for which the transition is a first order phase transition). They applied the usual 
periodic boundary condition (which approximates the thermodynamic limit in a very effective way) 
and also a disorder wall boundary condition. The difference between the T c values can be as large 
as 30 MeV (c.f. Figure [T^). Clearly, a detailed study of this phenomenon is needed for dynamical 
QCD. 
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5. The renormalized 

static potential in the continuum limit 

Our aim is to compute the free energy 
of a static quark-antiquark pair. There are 
several measurements on this quantity in the 



literature (for recent publications cf. [12, 13 



14, 15]). Here we go beyond these compu- 
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Figure 10: Transition temperature for the pure SU(3) 
theory as a function of the spatial size for Nt=4 and 
6. The lower, almost horizontal lines represent the re- 
sults with periodic boundary condition, whereas the up- 
per curves show the results with disorder wall boundary 
condition (Figure from Ref. [pT|). 



tations (for details see [16]), we use physical 
quark masses and perform a careful contin- 
uum limit extrapolation with the necessary 
renormalization procedure. 

The quark-antiquark free energy can be 
expressed as correlators of Polyakov loops: 



— F- 

e qq 



^<TrP(x)Tr/ 5t (x + r) 



(5.1) 

where r is a vector in the spatial direction, 
: 1 / (N t a) is the temperature and x runs over all the spatial lattice sites. P is the Polyakov loop. 
In full dynamical QCD at large distances it is favorable to generate a quark-antiquark pair from 



the vacuum, which then screens the color field between the two Polyakov loops []17[]. From this 
point (the string breaking scale) the lowest energy level will be insensitive of the position of the 
heavy quarks, resulting in a constant free energy. 

At finite temperature the above picture persists, but we can also have general expectations 
about the temperature dependence. Physically we expect that in a thermal vacuum it is easier to 
generate a quark-antiquark pair than at T = 0, since there are thermally excited particles around 
which can scatter on the gluonic string between the static quark-antiquark pair. The gluonic string 
can more easily break into a dynamical quark-antiquark pair. 

When we approach the continuum limit, the value of the un-renormalized free energy diverges. 
This is because in a single Polyakov loop the self-energy is divergent. We expect: 



(TrP(*)) 



e -C(a)N,a = e -C(a)/T 



(5.2) 



div 
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where C(a) —> °° in the continuum limit. At finite a the specific value of C(a) has no physical 
meaning, since it depends on how we define the "divergent part" of the self-energy (renormalization 
scheme). In the literature there are several ways to fix this constant [|l4|, ^]. A possible way of 
fixing C(a) is to take a physical observable based on Fq q , and require that it should be independent 
of a. In fact, the most useful quantity in our calculation was the constant value of the free energy 
after the string breaking/screening, at a fixed temperature (its value was set to Tq = 190 MeV 
in the calculation). The value of To = 190 MeV was motivated by the fact that it is already in 
the deconfined phase where the statistical errors of the free energy are much smaller than in the 
confined phase. For the calculation we used the gauge configurations of Ref [^]. Having determined 
the renormalized free energy for all lattice spacings, we could take the continuum limit by using 
the N t =4,6,8 and 10 free energies, and extrapolating in l/N 2 ~ a 2 —> 0. In the studied region 
N t = 8 and 10 results almost completely coincide. Therefore a safe extrapolation to l/N 2 = is 
possible. We estimate the systematic error of this extrapolation by comparing the results coming 
from N t = 6, 8, 10 extrapolation and N t = 8, 10 extrapolation. The result for the renormalized free 
energy at different temperatures, including both the statistical and the systematic errors, can be 
seen on Figure 1 1 . 
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of state at high temperatures 



While available lattice results for the 
equation of state (both for pure gauge 
theory and for full QCD) end at around 
5 • T c , standard perturbation theory con- 
verges only at extremely high tempera- 
tures (at 5 • T c , the different perturbative 
orders can not even tell the sign of the 
deviation from the Stefan-Boltzmann li- 
12 3 4 m jt). Until recently no link between the 

r [fm] 

two most systematic methods of QCD, 

r . „ rp, v ac ■ • t , namely perturbation theory and lattice 

figure 11: The renormalized tree energies in the continuum J v J 

v im j t QCD, existed for bulk thermodynamical 

quantities. 

It is easy to understand the difficulties for lattice calculations beyond these (e.g. 5 • T c ) tem- 
peratures. The correctly renormalized pressure is obtained by subtracting the vacuum term, thus 
Pren(T)=p(T)-p(T = 0), where p=T /V -log(Z) and Z is the partition function. Since the parti- 
tion function itself is difficult to determine, usually the integral method is used, which e.g. for the 
plaquette action is given by 

Pr,n(T)= P (T)-p{T = 0)= |<//3 «Pl) r - (Pl)o) (6.1) 

As an example imagine a calculation at 20 • T c on N t =% lattices. A back of an envelope estimate 
can convince you that it corresponds to a lattice spacings of about 0.0075 fm or lattice sizes of about 
1000 in full QCD. This is clearly out of reach in the near future. 
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In the perturbative approach the infrared properties of the non-abelian gauge theories at non- 
vanishing temperatures are responsible for the bad convergence. 

To create a link between these two system- 



ex 
a 



■ Boyd et al Nt=8 

O new approach Nt-6. 



atic approaches, we suggested a method [ ]18[ ] to 
obtain lattice results for the pressure at temper- 
atures, which were previously unreachable. It 
becomes possible to compare the data with per- 
turbation theory formulae. The results are ob- 
tained using two different techniques: the first 
one is a new way to renormalize the pressure, 
the second one is a direct method to measure 
the pressure. For illustration, the techniques and 
the results for pure SU(3) gauge theory are pre- 
sented. The extension to full QCD is straight- 
forward. 

The basic idea of the first technique is the 
observation, that the the divergences that are re- 
moved by renormalization are independent of 
the temperature. Thus the subtraction of the 
vacuum term, which is needed for renormaliza- 
tion, can be done at non- vanishing temperatures 
(actually even at very high temperatures). Thus, 
one can build up p mn as a sum of differences 
p ie n(THp(T)-p(T/2)]+[p(T/2)-p(T/4)]+. . . 
In fact, we usually measure the dimensionless pressure, which can be obtained by including the T 
factors. This leads to increasing powers of 1/2 4 in the successive terms, thus only a few of them are 
needed. Applying this scheme, we may reach arbitrarily high temperatures using lattices with only 
N t and 2 -N t temporal extents. It is worth mentioning, that a similar formula can be constructed for 
the case of the normalized interaction measure / = (e — 3 ■ p)/T 4 . 

Note, that independently of the renormalization procedure, a problem emerges within the in- 
tegral method framework. Since strictly speaking, the pressure is only exactly zero at T = 0, in 
principle we would have to carry out the integration starting from zero temperature. Thus, there 
is an uncertainty of setting the lower point of the integral. This problem is solved by our direct 
method (c.f. ^\). 

Let us consider the first bracketed term in p ien (T). 

\p(T)-p(T/2)\ = _i 7T logZ(^) - _l^iogZ(2tf,) = ^^log . (6.2) 



T/T 



Figure 12: The pressure, normalized to its Stefan- 
Boltzmann value, as a function of the temperature 
obtained by our new technique. Results with small 
discretization errors (Nt — 8, blue circles), seem to 
fit to the improved perturbation theory prediction, 
and also reproduce results obtained by the standard 
method at lower temperatures. At the highest tem- 
perature, 3 • 10 7 • T c , the pressure (within its statistical 
uncertainty) is consistent with the Stefan-Boltzmann 
limit. 



7 log , 

2N t N* \Z(2N,) 



N,Nf b y ' 2N,Nj & v ' 

This can be determined by using S\t (action with one boundary in the temporal direction), S%b 
(action with two boundaries in the temporal direction) and an interpolating partition function 
Z(a)=J 3)U exp[—(a-S2b + (l — cc)-Sib)] Using Z( a), one obtains 



W)-p(T/2)} 



lOE 



2(0) J 



1 , <flogZ(a) 
da — 
o da 



da{Si b -S 2b ) (6.3) 
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So we can calculate the pressure itself at any given temperature, without carrying out simulations 
at lower temperatures, and then performing an integral. (Note, however, that huge cancellations 
appear also within this framework). 

The method discussed above gives us the possibility to measure the pressure at very high 
temperatures. This was carried out using lattices with temporal extension N t = 4, N t = 6 and 
N t = 8. The temperature interval ranged from 4 • T c upto 3 • 10 7 ■ T c . Our results are shown on figure 



12 , A comparison can be done with the standard method, namely the result of []19(] for smaller 



temperatures. Here we present the results which are closest to the continuum limit, thus N t = 8. 
These results nicely follow the perturbative predictions ( for a recent paper see [ 2C ] and references 
therein). 



7. Conclusion 

QCD thermodynamics results from the Budapest- Wuppertal group were summarized. The 
necessary balance between T = and T > simulations was discussed in detail. As a consequence, 
Symanzik improved gauge and stout-link improved staggered fermionic lattice action was used in 
the simulations with an exact simulation algorithm. Physical masses were taken both for the light 
quarks and for the strange quark. The parameters were tuned with a quite high precision, thus 
at all lattice spacings the mx/fa and m^/m^ ratios were set to their experimental values with an 
accuracy better than 2%. Four sets of lattice spacings on lattices with _/V f =4,6,8 and 10 temporal 
extensions were used (they correspond to lattice spacing a~0.3, 0.2, 0.15 and 0.12 fm) to carry 
out the continuum extrapolation. It turned out that only A^=6,8 and 10 can be used for a controlled 
exuapolation, N t =4 is out of the scaling region. 

The nature of the T>0 transition was determined. The renormalized chiral susceptibility was 
extrapolated to vanishing lattice spacing for three physical volumes, the smallest and largest of 
which differ by a factor of five. This ensures that a true transition should result in a dramatic in- 
crease of the susceptibilities. No such behaviour is observed: the finite-size scaling analysis showed 
that the finite-temperature QCD transition in the hot early Universe was not a real phase transition, 
but an analytic crossover (involving a rapid change, as opposed to a jump, as the temperature var- 
ied). As such, it will be difficult to find experimental evidence of this transition from astronomical 
observations. Since for present day heavy ion experiments the baryonic chemical potential is also 
very small, the above results apply for them, too. 

The absolute scale for the T>0 transition was calculated. Since the QCD transition is a non- 
singular cross-over there is no unique T c . This well-known phenomenon was illustrated on the 
water-vapor phase diagram. Different observables led to different numerical T c values in the con- 
tinuum and thermodynamic limit also in QCD. Three observables were used to determine the 
corresponding transition temperatures. The peak of the renormalized chiral susceptibility pre- 
dicted r c =151(3)(3) MeV, whereas T c based on the strange quark number susceptibility resulted 
in 24(4) MeV larger value. Another quantity, which is related to the deconfining phase transition in 
the large quark mass limit is the Polyakov loop. Its behavior predicted a 25(4) MeV larger transition 
temperature, than that of the chiral susceptibility. Another consequence of the cross-over are the 
non-vanishing widths of the peaks even in the thermodynamic limit, which were also determined. 
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For the chiral susceptibility, strange quark number susceptibility and Polyakov-loop we obtained 
widths of 28(5)(1) MeV, 42(4)(1) MeV and 38(5)(1) MeV, respectively. 

The temperature dependent static potential was given. The same action and the same set of 
configurations were used as for the determination of T c . Since results for different lattice spacings 
were available a careful renormalization program was carried out. 

These features, numbers and functions are attempted to be the full result for the T^O transition, 
though other lattice fermion formulations — e.g. Wilson fermions (for ongoing projects see e.g. 
JH) , [22"| ]) or chiral fermions (for an early dynamical overlap test see [23], for the domain wall 
approach a recent presentation can be found in Ref.[24]) — are needed to cross-check the findings 
with staggered fermions. 

A new technique is presented, which -in contrast to earlier methods- enables one to determine 
the equation of state at very large temperatures. The method is based on the observation, that the 
divergences, which are removed by renormalization, are independent of the temperature. Thus the 
subtraction of the vacuum term, which is needed for renormalization, can be done at non- vanishing 
temperatures (actually even at very high temperatures). A direct method was also suggested, which 
does not need an integration over the temperature. Results for the pure SU(3) theory are presented 
upto 3 • 10 7 • T c temperatures. 
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